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ABSTRACT 

A series of 353 red electronic spectra (from three observatories, mostly from 6300 to 6700 A) obtained between 1994 and 2010, and 
of 171 UBV photometric observations (from two observatories) of the 2010 eclipse, were analyzed in an effort to better understand 
e Aur, the well-known, but still enigmatic eclipsing binary with the longest known orbital period (~27 yrs). The main results follow. 
(1) We attempted to recover a spectrum of the companion by disentangling the observed spectra of the 6 Aur binary failed, but we 
were able to disentangle the spectrum of telluric lines and obtain a mean spectrum of the F-type primary star. The latter was then 
compared to a grid of synthetic spectra for a number of plausible values of T^ff and log g, but a reasonably good match was not found. 
However, we conclude that the observed spectrum is that of a low gravity star. (2) We examined changes in the complex Ha line 
profiles over the past 16 years, with particular emphasis on the 2009-2011 eclipse period, by subtracting a mean out-of-eclipse Ho- 
profile (appropriately shifted in radial velocity) from the observed spectra. We find that the dark disk around the unseen companion 
has an extended 'atmosphere' that manifests itself via blueshifted and redshifted Ha 'shell' absorptions seen projected against the F 
star. Significantly, the Ho- shell line first appeared fiiree years before first contact of the optical eclipse when the system was not far 
past maximum separation. (3) Analyses of radial velocities and central intensities of several strong, unblended spectral lines, as well as 
UBV photometry, demonstrated that these observables showed apparent multiperiodic variability during eclipse. The dominant period 
of 66'!21 was common to all the observables, but with different phase shifts between these variables. This result strongly supports 
our earlier suggestion that the photometric variability seen during eclipse is intrinsic to the F star, and therefore, the idea of a central 
brightening due to a hole in the disk should be abandoned. Although variability on similar timescales is also seen in the spectrum and 
in photometry out of eclipse, we were unable to find a coherent periodicity in these data. Nevertheless, these results appear to rule out 
regular stellar pulsations as the cause of this variability. 

Ksy words. Stars: variables: general - Stars: binaries : eclipsing - Stars: individual: e Aur 



^> ■ 1. Introduction eclipsing obje ct must be huge. Analyzing the primary eclipse, 

k>( •, » • ^ T.T^ n^^^A TTr. i^oi-x ■ , ■ > iHuangI (Il965h concluded that the eclipsing object was a cool 

>< Epsdon Aurigae (7 Aur HD 3 1964, HR 1605) is a bright star- ^.^^ secondary eclipse has ever been observed. An 

^ ■ V^L'Z'fl'''-'''' °^ tT^' ^^Tu^: object in the center of the disk is invisible, and various hypothe- 

- ■ +43 49 24 ) It IS a smgle-line echpsing binairmh the longest ^^^^^ .^^ ^^^^ ^^^^^.^ ^^^^ ^^^^.^^ ^.^^^ ^^^^.^ 

toown orbital period of 27.1 years (989a?3) dl^ndorff 1903, ^^^-^^ individual masses and radii. It is there- 

I Chadimaetal .' 2010V Even more notable than the orbital period ^^^^ conceivable that 6 Aur belongs to those relatively rare bina- 

is the long duration of the eclipse of almost two yeai's, which ^.^^ ^^.^^ ^ ^yr) for which the brighter component in the opti- 

in combination with the l arge binary separation means that an ^p^^^^^j ^^^.^^ j^^^ ^^^^.^^ ^^^.^ 

Send offprint requests to: P. Chadima, confusion, throughout this paper we use the terms primary and 

e-mail: pavel.chadima@gmail.com secondary to denote the visible F-type star and the unseen object 

* Based on spectra obtained at the Dominion Astrophysical (presumably hidden in the cool disk), respectively. 

Observatory, Ondrejov Observatory and Castanet-Tolosan Observatory ^^^^ models of the system have been considered, 

PMx^oL™r^'^ ^''^^^'^'^ ^'^ ""^ Observatory and Hopkins pften referred to a s the high-mass, and the low-mass model 

f^"^^^, ■ ^ . ^T^,, jGuinan & DeWar jr2002i and references therein). 
Tables 1 and 2 are available only in electronic form at the CDS 

via anonymous ftp t o cdarc ■ u-strasbg . fr (130 .79. 128.5) or via The high - mass m odel was proposed in a remarkable study by 

|http : //cdsweb . u-strasbg. f r/cgi-bin/qcat? J/A+A/| . iKuiper et al.l (Il937h and consists of an FOIa primary of ~36M( 
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and a 1200-1400 K cool secondary hidden in the opaque disk, 
with a mass of ~2A.5Mq. They also estimated the radii of 
190 and 2690 Rq. The separation of the binary components is 
then ~35 AU. They pointed out that the observed flat-bottomed 
eclipse can only be understood if the cool eclipsing body is 
an extended gaseou s shell. The model was later developed by 
ICarroll et alJ (*1991) who considered the FOIa primary of - ISM© 
and the secondary hidden in the opaque disk, with a mass of 
~13.7M0. The separation of the binary components in their 
model is ~27.6 AU. The disk surrounding the secondary is con- 
sidered to have a protoplanetary origin. 

The low-mass model assumes a bloated post-AGB object 
and a secondary hidden in an optically thick disk, which is 
assumed to originate from recent accretion of matter flowing 
from the F primary to secondary. This model was first formu- 
lated by Eggleton & Pringle ( 1985) who suggested the masses 
of 1.3Mpi for the F star and 5Mq for the secondary. Along with 
iLissauer & Backmaa {19S4), they could not rule out the possi- 
bility that the object hidden in the disk is itself a binary, thus 
making e Aur a triple system. 

Although these two models are very different ones, there 
were no firm observational constraints to prefer one over the 
other. Very regrettably, the dista nces th at follow from the 
Hipparcos parallax (362-2273 pc: iPenym an & ESA 1997) or 
from the revised reduction (355-4167 pc; van Leeuwen 2007a b) 
are so impr ecise that they are consistent with both models. 
iHoard et"all (2010) assembled new Spitzer Space Telescope 
IRAC observations of e Aur , and after combining them with 
archival ultraviolet (UV) to mid-infrared data, they obtained a 
spectral energy distribution (SED) from -0.1 fim to 100 jjm. 
They argue that this SED can be reproduced by a three- 
component model consisting of a 2.2+ O.9M0 F type post-AGB 
star (the primary) and a 5.9+ O.8M0 B5V star (the secondary) 
surrounded by a geometrically thick but partly transparent disk 
of gas and dust with an effective temperature of 550 K. However, 
we note that the far-UV spectrum of e Aur observed by the 
FUSE satellite is an emission line spectrum, probably produced 
by scattering of continuum photons from a hot star embedded 
in th e occulting disk of the secondary (Bennett et al. 2005; Ake 
I2OO6 ). Future models of the spectral energy distribution should 
take this into account. 

Some additional su pport for the h i gh-ma ss model comes 
from the recent work of ISadakane et al.l (1201 Ol) . who performed 
an LTE abundance analysis of e Au r using A T LAS models and 
applying NLTE corrections. Sadak ane et al.l (1201 Ol) compared 
e Aur to the reference star HD 81471, an A7 lab supergiant. 
They find that both stars have comparable Tetf and log g, but that 
e Aur had higher micro turbulent velocity. Overall abundances 
of the elements Mg, Si, S, Ca, and TI, and of Sc, Cr, and Fe 
are comparable and close to solar in both stars. The i-process 
elements Y, Zr, and Ba are slightly more abundant (by +0.25 
dex) in e Aur than in HD 81471, but Sr is anomalously low in 
both stars. Carbon and oxygen are underabundant, and N and 
Na overabundant, in both stars. This abundance pattern is typ- 
ical of late-type supergiants and is quite different from that of 
post- AGB stars. In particular, for the four examined post-AGB 
stars. ISadakane et al. (2010) find enhanced carbon abundances, 
with [C/Fe] > [N/Fe] for these stars (in contrast to the situa- 
tion for e Aur and HD 81471). They conclude that the observed 
abundances of e Aur are normal for high-mass supergiant stars. 

It would be impractical to summarize all of the published 
studies of e Aur. We restrict ourselves to those relevant to the 
findings presented in this p aper. Several investigators, start- 
ing with iKuiper et al.l (Il937l) . have noticed the presence of a 



'shell' spectrum during the eclipse, which manifests itself by 
the presence of additional absorption lines, obviously originat- 
ing in the optically thin outer atmosphere of the disk. This 
shell spectrum can be us ed to analyze of the pro perties of the 
disk around the secondarv. lHinkle & Simonl(ll987l) observed CO 
shel l lines, which are absent in the spectrum of the primary, 
and Lambert & Sawyer (1986") found that the K I lines at 7664 
and 7699 A strengthen noticeably at the beginning and during 
the eclipse. Measuring their radial velocity (RV) and assuming 
a Keplerian disk, they estimated probable masses of the com- 
ponents of e Aur as less than 3 M©, and 3-6 Mq, for the pri- 
mary and secondary, respectively, and argued in favor of the 
low-mass model. They also observed cyclic line-profile varia- 
tions on a timescale of about two months and noticed that two 
known pulsating F-G la supergiants, HR 8752 and p Cas, have 
much longer periods of about 1 year. They speculated that, if the 
two-month variations are indeed caused by photospheric pulsa- 
tions of the primary, this ma y lend addit i onal s upport to the con- 
clusion that its mass is low. ISaito et al.l (Il987l) analyzed several 
shell lines and attempted to derive some properties of the disk. 
Since the shell lines were blended with the lines of the primary, 
they corrected the eclipse spectra, dividing them by one refer- 
ence spectrum taken out-of-eclipse. We note, however, that since 
the e Aur spectra are intrinsically variable, it may not be optimal 
to use only one, accidentally chosen reference spectrum for such 
a correction. Indeed, Ferluga & Mangiacapra (1991) used sev- 
eral out-of-eclipse spectra to derive an 'average' reference spec- 
trum, which they used for the reconstruction of 15 pure shell 
spectra. They then measured the physical quantities characteriz- 
ing the properties of a number of shell lines originating in the 
disk. Backman et al. (1985) analyzed hydrogen Brackett lines 
during the 1982-1984 eclipse. They concluded that the emission 
and absorption velocities of the Ba and By lines have a time de- 
pendence more like what is expected for the primary than the 
secondary, but with a large uncertainty. 

Readers are also r eferred to a review paper by 
iGuinan & DeW^ (I200I where detailed pieces of infor- 
mation about this mysterious system, as well as important 
references, can be found. 

Recently, a long and extensive series of archival radial veloc- 
ity beginning in 1899, and photometric observations of eclipses 
since 1848, have been compiled and analyzed independently by 
two groups. The results were the new, more precise ephemerides 
arid^orbital solutions of Stefanik et al. (2010) and Chadima et al] 
(l2010i) . These are clearly supe rior to the previous spectroscopic 
orbital solution of Wrightl (119701) since they are based on much 
longer observational periods and incorporate more accurate RV 
measurements. Moreover, they simultaneously use photometric 
measurements to derive a more accurate orbita l period and mide- 
clipse epoch. The ephemeris of IStefanik et al.l (1201 Ol) is 

rprim.min. = HJD (2455413.8 ± 4.8) + (9896'?0 ± 1^6) x E, (!) 

and it redicts the current mideclipse somewhe re between 
July 31 and August 10, 2010. IChadima et all (1201 Ol) derived the 
ephemeris 

Tpnm.min. = HJD (2455402.8 ± 1 .0) + (9890'?26 ± 0'?62) x E, (2) 

which predicts the current mideclipse to occur on July 24-26, 
2010. This orbital so lution is used in this paper. Furthermore, 
Chadima et al. (2010) argue that the brightening observed during 
the middle of the 1982-1984 eclipse was not due to a central 
hole in the disk, but was a continuation of the intrinsic F star 
variations observed out of eclipse. 



2 



p. Chadima et al.: Spectral and photometric analysis of e Aurigae prior to and during the 2009-2011 eclipse 



1.2 
1.1 

1 - 
0.9 - 
0.8 
0.7 
0.6 - 
0.5 - 
0.4 




Fe II 6433 
Fe II 6417 



Si II 6371 
Si II 6347 



H alplna 
J L 



6300 



6400 



6500 
wavelengtln A 



6600 



6700 



6800 



Fig. 1. Sample spectrum of e Aur from OND, showing the strong, unblended lines of Si II and Fe II for which the radial velocity, 
central intensity, and equivalent width were measured. 



The current eclipse began August 2009 and is predicted to 
end in May 201 1. Naturally, this rare event has attracted the at- 
tention of many astronomers in an attempt to better understand 
the nature of this still mysterious system. To this end, the present 
collaborative effort has obtained a diverse variety of observa- 
tions using the best available instruments and observational tech- 
niques and modern analysis meth ods. An outstanding example 
of such work is the recent study bv lKloppenborg et al.l (12010) in 
which the 2-D image of e Aur was reconstructed from CHARA 
multi-aperture, multi-baseline interferometry. Preliminary re- 
sults obtained during the current eclipse show an image of 
the primary star being partly eclipsed by a disk-like compan- 
ion, thereby providing direct confirm ation of th e Huang (1965) 
model. Assuming a mass of 5.9 (Hoa rd et al.i i2010) for the 
B star hidden in the disk, they estimated the mass of the F star as 
3.6 ± 0.7 Mb . 

iLeadbeater & Stence il (I2010h report a step-wise increase in 
the line strength of the KI line at 7699 A during the current 
eclip se, which seems to support an earlier suggestion by Ferlugal 
( 11990) that the disk has a multiring structure. These a nd other 
recent findings about e Aur have been summarized by IStencell 
(120101) . 

The present study is organized as follows. The observational 
material and its reduction and measurement are presented in 
Sect. |2l Our attempts at disentangling the e Aur spectra and the 
comparison of disentangled spectra of the primary with the syn- 
thetic spectra are discussed in Sect. |3] Analyses and possible in- 
terpretation of the complex and highly variable Ho' line profiles 
are given in Sect.|4] An analysis of the measured radial velocities 
and central intensities of several spectral lines and UBV photom- 
etry is presented in Sect. |5] 



2. Observational material, reduction, and 
measurements 

2.1. Spectroscopy 

This study is based on the following three series of electronic 
spectra covering the spectral region around the Hot line. 

1. 105 CCD spectra with a linear dispersion of 10 A mm ' and 
a two-pixel resolution of 21700 were secured in the Coude 
focus of the 1.22-m reflector at the Dominion Astronomical 



Observatory (DAO) in Canada between 1994 and 2010 by 
S Y and PDB. These observations cover the wavelength range 
6200-6750 A. 

2. 201 CCD spectra having a linear dispersion of 17 A mm ' 
and resolution of 12600 were obtained in the Coude fo- 
cus of the 2.0-m reflector of the Ondfejov Observatory 
(OND) in the Czech Republic in 2006-2010 by PH, PS, 
MS, MW, PC, VV, and several additional observers (credited 
in Acknowledgements). These observations cover the wave- 
length range 6260-6760 A. 

3. 47 CCD spectra with resolution of 11000 were obtained 
with a 0.28-m Celestron 1 1 telescope at the Castanet-Tolosan 
Observatory in France (CTO hereafter) by CB during 2009- 
2010. 

Initial reductions of the DAO and OND spectra were carried 
by SY and MS, respectively, both in IRAF. The initial reductions 
of the CTO spectra were carried out by CB in Reshel V 1 . 1 1 . 
Rectification and all radial velocity, central intensity (CI), and 
equivalent width (EW) measurements of sele cted lines in these 
spectra were carried out by PC in SPEFO (iHorn et al.lll996l 
ISkodalll996h . The latest JK 2.63 version of die program SPEFO, 
developed by the late Mr J. Krpata, was used. All measurements 
are presented in detail, together with heliocentric Julian dates 
(HJDs) of mid-exposures in Table 1 (electronic only). 

A representative Ondfejov red spectrum is displayed in 
Fig- [H Quantitative measurements were carried out for sev- 
eral of the stronger, unblended spectral lines: two Si II lines 
at 6347.109 A and 6371.371 A, two Fell lines at 6416.919 A 
and 6432.680 A, and the Ho- line at 6562.817 A, for which the 
emission wings were also measured. The RV measurements in 
SPEFO are based on sliding the direct and reverse line profile on 
the computer screen until the best coincidence of both profiles is 
achieved. The zero point of the RV scale was determined indi- 
vidua lly for each spectr um via measurements of selected telluric 
Unes (lHometal.lll996h . 

The behavior of the Ha profile is illustrated well by the 2-D 
images of the line profile over time, shown in Fig. |2]for the entire 
1994 to 2011 period, and in Fig.|3]in detail for the 2009-2011 
eclipse. Note the high degree of variability of the Ha profile. The 
velocity scale for these spectra is heliocentric, and one can see 
the annual sinusoidal variations of several nearby telluric water 
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Fig. 2. Evolution of the Ha profile from 1994 to present. The x-axis shows heliocentric velocities (in km s ') and the y-axis shows 
observation dates (HJD-2,400,000). The intensity scale is shown below. The spectra were linearly interpolated onto a regular time 
grid to produce this image. The black dashes along the right axis show the actual observation dates. 
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Fig. 3. Evolution of the Ha profile during the 2009-201 1 eclipse. The x-axis shows heliocentric velocities (in km s ') and the y-axis 
shows observation dates (HJD-2,400,000). The intensity scale is shown below. The spectra were linearly interpolated onto a regular 
time grid to produce this image. The black dashes along the right axis show the actual observation dates. 
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vapor lines in Fig. |3] A detailed analysis in the Ha variations is 
presented in Sect. |4] 

2.2. Photometry 

Here, we use the following two series of photoelectric UBV ob- 
servations from the phases of the current total eclipse. 

1. 105 differential observations relative to A Aur (HD 34411) 
in a period HJD 55213-55468 were obtained by HB, DS, 
PH, MW, and PC with a photoelectric photometer and an un- 
cooled EMI 9789QB tube at the Hvar Observatory (Croatia) 
0.65-m Cassegrain reflector. Observations in all three filters 
consisted of from 7 to 10-second integrations followed by 
5-second integrations on the sky. The check star HR 1644 
(HD 32655) was observed as frequently as the variable and 
typically three to five such individual differential observa- 
tions were obtained for each night of observations. All these 
observations were reduced and transformed to the standard 
UBV system with the help of non-linear transformations us- 
ing the reduction prog ram HEC22 (Harmanec et al. 1994; 
iHarmanec & Horn 1998). The latest release #17 of the pro- 
gram HEC22 was used. This version allows modeling the 
varying extinction that occurs during the course of each 
night's observations. The standard UBV magnitudes of the 
comparison star A Aur 

V = 4"?705, B-V= 0'%21, and U-B= 0"?142 

were derived from numerous all-sky observations at Hvar 
and were added to the magnitude differences var-comp and 
check-comp to obtain final differential UBV magnitudes of 
e Aur and HR 1644. All Hvar observations of both stars were 
tabulated with the HJDs in Table 2 (electronic only). 

2. 66 differential observations, also obtained relative to A Aur 
in a period HJD 55200-55539 were secured by JH at the 
Hopkins Phoenix Observatory, USA. His home-built pho- 
tometer uses a 1 P2 1 photomultiplier tube operated at -950 V 
and is attached to a Celestron C-8 telescope mounted on a 
Meade LX-90 mount. Dead time for the detector is deter- 
mined using an aperture mask. The mask contains irregu- 
larly spaced holes that reduce the aperture by 80%. Stars 
that are close to each other and near the zenith are measured 
as quickly as possible to reduce the effects of extinction. A 
bright star and a faint star are measured without, with, and 
without the aperture mask placed in the light path. Their re- 
sulting counts are averaged and used to compute the dead 
time coefficient of the instrument. Observations are con- 
ducted in the standard CVCVCVC format (C = comparison, 

V = variable). Each observation is composed of three 10- 
second integrations in each of the three UBV filters, followed 
by one 10-second integration on the sky, again in UBV. The 
following values of the comparison star, slightly different 
from the Hvar values, were adopted: 

V = 4™71, B-V= 0™63, and U-B= 0"?12. 

These datiO are on the instrumental UBV system of the in- 
strument. 

3. Disentangling and modeling the e Aur spectra 

Having at our disposal several hundred e Aur spectra, which 
span an interval of 17 years, i.e. more than a half of the orbital 

' Available from: 

|http : //www.hpsbft . com/EAUrO 9/Data/UBVRIData . html| 



period, we attempted to disentangle them in the hope of detect- 
ing a faint spectrum of the 'unseen' secondary, possibly lost in 
the much more prominent spectrum of the primary. Note, for in- 
stance, that if the true mass ratio of the system was not too far 
from one, the lines of both components would remain blended in 
all orbital phases since the semi-amplitudes of their RV curves 
would both be only -15 km s ' . For the purposes of disentan- 
gling spectra, we only used observations obtained prior to the 
predicted beginning of the current eclipse. This left us with 208 
spectra (99 DAO and 109 OND spectrograms). We mainly used 
the interval 6333-6540 A and also several shorter subintervals 
(e.g. an interval of 6344-6374 A, containing both Si II lines). 
We made no attempt to disentangle Ho- , considering its compli- 
cated shape, obvious VI R changes of its emission wings, and the 
possibility that it need not be associated with either the primary 
or secondary component. 

We used the program KOREL, developed by Ha draval (I1995L 
1 19971 1200 4^ and made publicly available by the author. It uses a 
spectral disentangling technique in the Fourier domain and sim- 
plex method of minimalization of the sum of residuals and dis- 
crete spectra from various orbital phases, which are rescaled to 
a linear RV scale. From these data, the spectra of the individ- 
ual binary or multiple components are derived, as well as, if re- 
quired, the orbital elements of the binary or multiple system. No 
assumption is made about the shape of the line profiles; the only 
requirement inherent to the method is that the component pro- 
files can only vary in their strength and not in their shape. The 
run of KOREL is controlled by a parameter file with the initial 
orbital parameters of the system and initial simplex steps. Any of 
the orbital parameters can either be held fixed or be converged. 

We selected only a short wavelength interval 6529-6539 A 
containing virtually no stellar lines but only a number of stronger 
telluric lines, then fixed the orbit of the Earth around the Sun pro- 
jected into the direction to e Aur, and let the program derive the 
relative strength of the telluric spectrum for each spectrogram 
used. Then we fixed these strengths and also all elements of the 
orbital solution derived by Chadima et al. (2010). We instructed 
the program to proceed step by step, first allowing for the deter- 
mination of variable line strengths for both binary components 
only and finally also allowing for convergence of the epoch of 
periastron, which is the most uncertain parameter in the adopted 
orbital solution. 

To see whether the resulting sum of squares of residuals 
is sensitive to the mass ratio q used (thus indicating the definite 
presence of a secondary line spectrum), we carried out a number 
of KOREL solutions for different values of q. Such a mapping 
has been successfully applied i n several other cases - see, for 
instance. iHarmanec et al.l(l2010h . 

We applied this mapping to the whole selected wavelength 
interval, as well as to several sub-intervals but the expected, 
roughly parabolic dependence of on the mass ratio has ever 
been found. We must conclude that, at least for the spectral reso- 
lution we had at our disposal, no lines of the secondary spectrum 
were detected. 

After obtaining this null result, we decided to use KOREL 
only for disentangling the telluric spectrum and recovering the 
mean out-of-eclipse spectrum of the F-type primary. Using the 
mean stellar spectrum minimizes the effect of transient, intrinsic 
variations. 

We used a catalog of model atmospheres ATLAS for our cal- 
culation of synthetic spectra. We calculated a grid of the syn- 
thetic spectra of 'classical' F stars based on Kurucz' ATLAS9 
model atmospheres for various temperatures, gravities, and 



6 



p. Chadima et al.: Spectral and photometric analysis of e Aurigae prior to and during the 2009-2011 eclipse 



1.2 
1.1 
1 

0.9 - 
0.8 - 
0.7 - 
0.6 
0.5 




mean, out-of-eclipse spectrum 
Tgff = 8000K, log(g) = 1 .0[cgs], log Z/X = 0.3 
Tgff = 8000K, log(g) = 3.0[cgs], log Z/X = 0.3 

T„ff = 8000K, log(g) = 1 .0[cgs], log Z/X = -0.3 

I I I I 



6340 



6360 



6380 6400 
wavelength A 



6420 



6440 



Fig. 4. Comparison of the mean, out-of-eclipse spectrum of e Aur with some ATLAS synthetic spectra. The best-fitting synthetic 
spectrum has the parameters: T^s = 8000 K, log g = 1.0[cgs] log Z/X = -0.3, v(rot) = 5.0 km s"'. Two other synthetic spectra are 
shown, one with higher gravity and one with lower metallicity. 



metallicities. In particular, the grid was derived for T^s = 7500- 
8500 K, log g = 1.0-5.0[cgs], and metallicities log Z/X from 
-0.3 to H-3.0. All synthetic spectra were convolved to rotation 
velocities of the primary in an interval v(rot) = 5.0-25.0 km s"' . 
Synthetic spectra for higher metallicities and lower gravities 
were not modeled because the ATLAS catalog does not con- 
tain model atmospheres for those parameters owing convergence 
problems. 

The best match between the disentangled and synthetic spec- 
trum (in the sense of the lowest value of x^) was obtained for 
a synthetic spectrum with the following characteristics: T^ff = 
8000 K, log g = LO[cgs], log Z/X = -0.3, v(rot) = 5.0 km s"'. 
Our results are in fair agreement with the spectral analysis by 
ISadakane et al.l (12010) . In Fig. 31 the disentangled spectrum is 
compared with several synthetic spectra, and it is clear that even 
the model for log g=l.Q does not describe the observed spectrum 
properly; that is to say, the two Si II lines (on the left side of the 
plot) and two Fe II lines (on the right side of the plot) should 
be more prominent than they are for the best synthetic spectrum 
with log g= 1 .0. 

To check on the result, we investigated the variation inx^ as a 
function of the input quantities, which characterize each model. 
We also inspected the plots of individual spectra. The conclu- 
sion from this exercise is that the synthetic spectra are almost 
insensitive to the value of the projected rotational velocity for 
the values used, at least in comparison to other characteristics. 
This implies that the 'optimal' value of the projected rotational 
velocity should not be given too much weight, even more so be- 
cause a satisfactory fit between the disantangled and model spec- 
tra was not found. The dependence on T^^fi is almost negligible 
for the range of investigated values and the spectral lines used. 
There is, however, a very strong dependence on the gravity and 
metallicity. In Fig.lH there are a few synthetic spectra for lower 
metallicity and higher gravity. It is seen that one could probably 
reproduce the Si II and the Fe II lines by increasing the metallic- 
ity and further decreasing log g . Regrettably, we do not have a 
suitable model atmosphere available to confirm this conjecture. 

However, the best-fitting synthetic spectrum already pro- 
vides a good fit to the weaker lines from 6375-5415 A, even 
though it fails to adequately describe the stronger Si II and Fe II 
lines. This behavior illustrates the limitations of the ATLAS 
model atmospheres used here, because they are plane-parallel 



LTE models. Both of these assumptions are probably signifi- 
cantly violated in low-gravity, highly evolved stars (regardless 
of whether the F star is of low or high mass). Therefore, these 
ATLAS models should be considered as providing only qualita- 
tive constraints on the stellar parameters (Teff , log g , abun- 
dance). For quantitative work, non-LTE (NLTE) models, and 
probably spherically symmetric models also, are needed. 

4. Analysis of the Ha profile during the 
2009-201 1 eclipse 

The profile of the Ha line of e Aur observed outside the eclipses 
consists of a central absorption surrounded by a double emission 
peaks, which might or iginate in the extended envelope of the pri- 
mary (ICha et al.ll 19941) . The whole Wa profile and, especially, its 
emission wings are very variable, which is sometimes explained 
by the rotation and physical changes in the envelope around the 
primary. This variability causes problems in the interpretation of 
the i\a profile even out of eclipse. 

To investigate the assumption that the Ha emission origi- 
nates in a vicinity of the primary, we calculated orbital solutions 
for the Si II lines and the Ha emission wings using RV measure- 
ments on the DAO out- of-eclipse spectra. We used the FOTEL 
pr ogram (iHadrava 2004) with input parameters from the solution 
of IChadima et al.l (12010) and only allowed a convergency of the 
periastron epoch T and the velocity semi-amplitude K\ . Results 
are summarized in Table |3] 

In Fig.|5] the orbital solution found from the Si II lines alone 
and the solution from the wings of Ha are compared with the 
solution of Chadima et al. (2010). It is seen that the Ha emis- 
sion roughly follows the orbital motion of the Si II lines, i.e., the 
primary, and therefore must be associated with this component. 
We infer that a small difference between both solutions is mainly 
caused by an insufficient interval covered by data used for both 
orbital solutions (when compared with the orbital period of the 
e Aur). 

To verify that the Ha emission originates near the primary, 
we selected two rather symmetric Ha profiles from substantially 
different orbital phases and overplotted them on a heliocentric 
wavelength scale in Fig. |6] One can see that the whole Ha pro- 
file, including the double emission wings, undergoes the RV shift 
corresponding to the motion of the F star in orbit. This behavior 
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Table 3. Orbital solutions for T, K\ , vq from fits to the Si II lines 
and Ho- wings. 
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Fig. 5. Orbital solutions derived only from the Si II lines and 
from the emission wings of Ha are compared to the full 
IChadima et"!!] (l2010h solution. 
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Fig. 6. Comparison of two symmetric Ha profiles from different 
orbital phases. The orbital phases and the radial velocities of the 
primary at those phases are labeled. 
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Fig. 7. Variation of the Ha core central intensity in the DAO and 
OND spectra. 



provides strong evidence that the Ha emission is associated with 
the F-type primary. 

During eclipse. Ha undergoes prominent line profile 
changes. Figure |7] shows a remarkable fact: the CI of the Ha 
absorption core steadily decreases with the approaching eclipse. 
An important finding is that this decrease started about three 
years before the beginning of the optical eclipse, thus confirm- 
ing the behavior first reported by iKuiper et al.l (Il937l) . The Ha 
EW increased during the same time period. The RV of the cen- 
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tral absorption core increased significantly during the first half 
of the eclipse, but after mideclipse, then suddenly reverted back 
to negative values. The Ha emission peaks disappear during 
eclipse. For the current eclipse, the red emission peak had dis- 
appeared by HJD 2455259 and the blue peak had vanished by 
HJD 2455329. The red emission peak subsequently reappeared 
after HJD 2455524. 

From a study of the Ha line profile during the last eclipse 

(iBar sonv et al. 1986), we expect the evolution of the Ha profile 
in the second half of the current eclipse to follow the sequence of 
events described above in reverse. The observations should show 
a gradual increase in the Ha absorption core CI, accompanied by 
a decrease in its EW, and a gradual return of the absorption core 
RV to the out-of-eclipse values. It is expected that both emission 
peaks will reappear in reverse order to what was seen in the first 
half of the eclipse. 

Based on the results of previous st udies of the shell spe ctrum 
(discussed in Sect.[T]l, we adopted the Kuiperet alJ(ll937h inter- 
pretation of the Ha behavior We assume that (apart from any 
intrinsic line profile changes, as seen for other lines) changes in 
the Ha profile during eclipse are due to additional absorption 
of the primary's light by the outer, optically thin atmosphere of 
the companion's occulting disk. In other words, we suppose that, 
since the companion's disk does not contribute any radiation in 
the optical region near Ha, the uneclipsed part of the primary 
star produces a similar Ha profile to the whole stellar disk out 
of eclipse. However, the region producing the Ha emission may 
be extended, in which case much of the emission may not be 
eclipsed during the stellar eclipse, thereby also resulting in a 
similar eclipse emission profile to the one out of eclipse. We in- 
terpret the disappearance of the out-of-eclipse Ha emission dur- 
ing eclipse (Fig.|2]i as the result of strong shell absorption at the 
corresponding wavelength. In this scenario, the intrinsic emis- 
sion is still present, but is strongly absorbed by material in the 
occulting disk during eclipse. This disk absorption reflects the 
velocity field in the rotating disk, which is why its RV increases 
in the first half of the eclipse. This naturally leads to an apparent 
disappearance of the red emission peak. As the eclipse goes on, 
the shell absorption from the other edge of the disk, rotating in 
the direction towards the observer, becomes dominant and the 
observed RV of the deepest absorption quickly shifts to negative 
(i.e., blueshifted) values, thus leading to an apparent disappear- 
ance of the blue emission peak. Near the center of the eclipse, 
both edges of the rotating disk atmosphere are projected against 
the primary and one should simultaneously see a blueshifted and 
a redshifted absorption. If the disk were axially symmetric, one 
should see two absorption cores of the same CI with opposite 
RVs with respect to the instantaneous RV of the secondary. 

To test and explore our hypothesis, we once again used 
KOREL to obtain a mean Ha profile, free of telluric lines and in- 
trinsic variations, using only the spectra observed out of eclipse. 
This mean Ha profile was then shifted to the primary's RV, using 
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Fig. 8. Evolution of the observed Ha profile near mideclipse. 



the IChadima et alj (1201 Ol) orbital solution, and subtracted from 
the observed OND spectra. Since the shell absorption should 
have the out-of-eclipse Ha as its continuum, we finally renor- 
malized the subtracted spectra to a continuum value of unity. In 
Fig. [8] we show observed Ha line profiles at five epochs near the 
date of mideclipse. The complex Ha line profiles shown here 
are similar to those reported by other observers during the cur- 
rent e Aur observing campaign. In Fig. |9] we show the mean, 
out-of-eclipse profile and the same five profiles as in Fig. |8] but 
with the mean, out-of-eclipse profile subtracted from them. The 
result of the subtraction supports our hypothesis and the valid- 
ity of our procedure. This is because the subtracted spectra show 
simple, double-absorption core profiles consistent with what is 
expected from the occulting disk model. Furthermore, these Ha 
absoiption profiles evolve during the course of the eclipse in a 
manner consistent with what is expected from the occulting disk 
model, in contrast to the complex observed Ha profiles of Fig. |8] 
The red absorption core is much deeper on the earliest spectrum, 
when most of the light from the primary is absorbed by a reced- 
ing part of the rotating disk. Furthermore, the ratio of the central 
intensities of the red/blue absorption cores gradually decreases 
through unity near mideclipse. At mideclipse (the third spectrum 
in both figures), the two absorption cores are almost symmetric 
(after subtraction of the mean, out-of-eclipse profile). 

Similar profile changes should also be seen for the li/3 line. 
We checked 30 available OND spectra covering the li/3 region, 
observed between 2007 and 2010. Regrettably, we have no out- 
of-eclipse H/3 spectra and cannot repeat the procedure carried 
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Fig. 9. Evolution of the corrected Ha profile, with the mean out- 
of-eclipse profile subtracted, near mideclipse. 



out for Ha. We therefore only present a selection of five original 
spectra in Fig.[TO] It is seen that their CI is decreasing and the ab- 
sorption core is shifting redward as the mideclipse approaches. 
Duiing and after the mideclipse, the absorption core quickly be- 
comes blueshifted. This indicates that the evolution of the H/3 
profile near the current mideclipse is indeed similar to that of 
the Ha profile. 

We measured the RVs and CIs of both Ha cores after the sub- 
traction of the mean, out-of-eclipse profile, as well as the EWs 
of the entire double absorption profile. All these quantities are 
plotted in Fig.[TT] It is seen that the CI decreases and the EW in- 
creases with the approach of mideclipse. These changes become 
more prominent closer to mideclipse as the light from the pii- 
mary goes through denser parts of the occulting disk, resulting 
in stronger disk absorption. The RV increases as the line-of-sight 
absorption through the occulting disk gradually moves closer to 
the disk's axis of rotation (and to correspondingly higher orbital 
velocities). The blue core, with its RVs in the opposite sense, first 
appears 74 days before the mideclipseQ, and it deepens quickly, 
while the the red core reaches its maximum strength and gradu- 
ally weakens thereafter The red core finally disappears 49 days 
after mideclipse, at about the time when the blue core attains 
its maximum strength (i.e., minimum CI). The interval when the 
both cores are visible is therefore 123 days. This is in a good 



- The difference between the HJDs of mideclipse and the first obvi- 
ous appearance of the blue core in the spectrum. 
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Fig. 10. Time evolution of the original Hf) profile during the 
2009-2011 eclipse. 



agreement with the estimated time needed for the center of the 
companion's occulting disk to transit the disk of the F star It 
is also seen that, while the EW of the entire Ha absorption line 
profile reaches its maximum several days prior to mideclipse, the 
red and blue absorption cores do not attain equal CI until several 
days after mideclipse. One possible explanation of this finding is 
that the center of the disk is not perfectly axially symmetric and 
that there is a denser region in the receding part of it and a more 
rarified region in its approaching part. It is important to continue 
these Ha observations for the second half of the current eclipse 
to confirm the occulting disk model proposed here. 

We emphasize that the RVs in Fig. [TT| should be corrected 
for orbital motion before any analysis. But since the 'shell' ab- 
sorption cores originate in the disk around the secondary, their 
RV measurements should be corrected for an orbital velocity of 
the secondary, not the primary. Nevertheless, uncertainties in 
present estimates for the mass ratio are too high, thereby keeping 
this correction from being inconvincing. However, it was proven 
that any reasonable value of the mass radio does not change 
the resulting RVs so much that it would affect our conclusions 
presented in this section. However, some inconsistencies remain 
with the occulting disk scenario. 

1. The average between the first (HJD 2455329) and 
the last (HJD 2455452) corrected spectrum for which 
both cores are visible is HJD 2455391, which dis- 
agrees with the predicted times of the mideclipse of 
either HJD 2455413.8+ 4.8 ( Stefanik etafl l2010l) or 
HJD 2455402.8+ 1.0 (C hadima et aL201Q) . 

2. The RVs of the blue core are gradually decreasing during last 
few observations, while their increase is expected. 

3. The last few Cls of the blue core (of the subtracted spectra) 
are slightly negative, which is, of course, physically impos- 
sible. 
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Fig. 11. Changes of RV (top), CI (middle), and EW (bottom) of 
the two absorption cores of Ha during eclipse, after subtraction 
of the mean, out-of-eclipse Ha profile. The first vertical line de- 
notes the beginning of the current photometric eclipse, the sec- 
ond one indicates the date of mideclipse. 



All of these inconsistencies can be explained if we realize that 
the out-of-eclipse Ha profile is highly variable, and it is un- 
doubtedly also variable during the eclipse, thereby affecting 
the observed eclipse profiles. Since we corrected all Ha pro- 
files by subtracting the mean profile, any intiinsic line variations 
will still affect the individual profiles Q Moreover, the measured 
quantities in the absorption profiles could also be partly affected 
by relatively strong telluric lines at 6560.555 A and 6564.206 A. 

Our hypothesis and results discussed above are potentially 
good for determining the physical properties of the disk (i.e., 
its density and velocity structure), which could help in under- 
standing its nature. But for a precise analysis of this kind, one 
should model the light of the primary going through the vari- 
ous model disks and compare the results with the observed Ha 
profiles. This task is beyond the scope of this, rather qualitative, 
analysis. 

4. 1. Sightlines at the onset of spectroscopic and 
photometric eclipse 

As already noted at the beginning of this section, the CI of the 
Ha absorption core starts to decrease at least three years be- 
fore the predicted beginning of primary eclipse. This observa- 
tional result can be used to constrain the structure of the system. 
We refer to the peiiod when Ha absorption cores deepen be- 
low their out-of-eclipse mean value as 'spectroscopic eclipse'. 
From Fig.|7] we estimate that spectroscopic eclipse began around 
HJD 2454000. The Chadima etal. (2010) orbital solution im- 
plies that the angle between the sightline at the start of spec- 



^ Note that lCha et al.l ( 119941) made an unsuccessful attempt to find pe- 
riodic behavior in the Hor profile; see also the first paragraph of Sect.|5] 
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Fig. 12. Critical Roche lobes for the high and low-mass models 
are shown to scale, with the lines of sight to Earth shown for the 
beginning of spectroscopic and photometric eclipse. 



troscopic eclipse and the mideclipse sightline is ~85°. JH deter- 
mined that the 2009 photometric eclipse began on HJD 2455056 
when the sightline angle was only ~20° relative to the mideclipse 
line of sight. 

In Fig. [12] we show the critical Roche lobes for the high- 
and low-mass models presented in Sect. [T] and the sightlines at 
the onset of spectroscopic and photometric eclipse. It is evident 
that the additional Hot absorption at the start of spectroscopic 
eclipse cannot be caused by material near the secondary because 
the binary orientation is close to maximum separation at that 
time. Therefore there must be some circumbinary mate rial re- 
sponsi ble for this additional absorption, suggestive of the lStruv3 
model. However, from our observations, it is seen tha t 
this envelope is not homogenous, as proposed bv lStruve1(ll956h . 

The line of sight at the onset of photometric eclipse inter- 
sects both critical Roche lobes around the secondary. However, 
in the case of the high-mass model, the disk around the sec- 
ondary would almost extend to the critical Roche lobe. 

4.2. A remarkable out-of-eclipse Ha profile variation in 
2005-2006 

An unusual and prominent change in the Ha line profile was 
observed during 2005-2006, long before the onset of spectro- 
scopic eclipse. The blue emission wing disappeared and was 
replaced by a deep, blueshifted absorption core (Fig. [T3] also 
Fig.|2]l. Initially, the first spectrum (in Apr 2005) has a 'standard' 
Ha profile, while the second and the third spectra observed later 
that year show the gradual appearance of the blueshifted absorp- 
tion core. By Mar 2007, the final spectrum of this series shows 
a return to the standard Ha profile, but now with an increase 
in EW due to the onset of spectroscopic eclipse. Our observa- 
tional coverage of this period is not sampled well enough to per- 
mit a reliable estimate of the timescale of the event. However, 
considering the first and the last spectrum for which a devia- 
tion from the standard Ha profile was seen, we estimate that the 
event lasted from August 2005 to March 2006. A similar behav- 
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Fig. 13. Dramatic Ha profile variation: the sudden appearance 
of a deep, blueshifted absorption core, as observed at the DAO 
during 2005-2006. 

iorof the Ha profile was also observed and reported bv IS channel 
(l2007l) . However, a similar profile variation is not observed for 
other spectral lines. 

The unusual change in the Ha profile was most probably 
caused by additional strong absorption from gaseous material 
with at an RV of about -40 km s"' relative to the primary. This 
absorption feature indicates the presence of either a transient 
outflow of material from the primary star or of some localized 
circumstellar matter producing an 'eclipse' of the primary. 

5. A period analysis of thie spectroscopic 
quantities and UBV phiotometry during thie 
2009-2011 eciipse 

The profiles of the selected spectral lines, measured by us, are 
all variable on much shorter timescales than for the binary or- 
bital period. To find out whether this variability is irregular or 
(multi)periodic, we carried out period analyses for individual ob- 
servables, using two different computer programs. The first pro- 
gram was the HEC2fl a code developed by R Harmanec an d 
based on phase dispersion minimalization ( Stelling werfll 1 978h . 
The second program, PERIOD0 4, developed by R Lenz and 
M. Breger (iLenz & Bregerl 120051) . uses Fourier analysis meth- 
ods; i.e., the time series of observational data is fitted using a 
sum of cosine functions. 

Initially, we analyzed the periodic behavior of the Ha profile 
out of eclipse. For analysis of the Ha absorption profile RVs and 
CIs, we selected the 84 DAO observations prior to HJD 24549 1 3 . 
Later observations were excluded because these were affected 
by the peculiar line profile behavior discussed in Sect. I4.2l and. 
subsequently, by the onset of spectroscopic ecUpse discussed in 

The program with brief instructions how to use it is available at 
http ://astro.troja. mf f . cuni . cz/ f tp/hec/HEC2 7, . 
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Fig. 14. RVs of the Fe II lines (upper panel), the CIs of the 
6416.919 A Fe II line (middle panel) and the B photometric filter 
(lower panel) fitted with two sinusoids. In each case the dom- 
inant period is 66.21 days, while the secondary period differs 
among the three cases. 



Sect. 14. II All out-of-eclipse measurements were included for an- 
alyzing of the emission wings: the RVs, the CIs of the red and 
the blue emission, the emission-peak ratio V/R, and the emission 
strength {V + R)/2 where V is the peak intensity of the blue wing 
and R is the peak intensity of the red wing. Both the HEC27 and 
PERJOD04 analyses found no evidence of periodic behavior in 
an y of the observab les analyzed, a result in accordance with that 
of ICha et al.l (Il994l) . We conclude that the physical changes in 
the circumstellar matter responsible for the Ha emission are not 
periodic. 

Next, we analyzed the behavior of the RVs and the CIs of 
the two measured Fe II and Si II lines. Since we assumed that the 
ecUpse could affect these lines and since we did not have enough 
observations during the eclipse, we initially analyzed only the 
out-of-eclipse data. We used both methods but we also did not 
find any periodic behavior in the changes of these lines. 

After collecting sufficient observations of the current (2009- 
2011) eclipse, we repeated the periodic analyses just for the 
eclipse period, motivated by indications of periodic behavior 
during eclipse. We used the PERIOD04 code for this analysis, 
and we found that all the observables could be reasonably mod- 
eled by two sinusoids, except for the RVs of the Si II lines, which 



e Aurigae prior to and during the 2009-2011 eclipse 

Table 4. Periods found during the 2009-201 1 e Aur eclipse. 



Observable 


Period 


Phase shift 




(d) 




Fe II and Si II RVs 


66.21 


0.000 




123.3 


0.668 


(just for Si II ) 




U.io/ 


Fe II CIs 


66.21 


0.808 




119.9 


0.659 


Si 11 CIs 


66.21 


0.808 




268.2 


0.356 


UBV photometry 


66.21 


0.323 




269.8 


0.081 



needed three harmonic functions to adequately describe their 
variations. It is interesting to note that the main period of ~66 
days was common to all the analyzed observables. In contrast, 
the secondary period differed among the observables. However, 
subsets of the observables shared common secondary periods: all 
the RVs had the same secondary periods, as did each of the Fe II 
line CIs, and each of the Si II line CIs (Tabled. As a final ex- 
ercise, we tried to extrapolate these periods to the out-of-eclipse 
observations (which have much sparser time coverage), but this 
failed to accurately describe the observed variability. There is 
evidence that the photometric variability of e Aur is quite coher- 
ent at some times, but this coherence is lost on longer timescales. 
This 'unpredictable' behavior is presumably due to some intrin- 
sic variability of the F-type primary star. 

We decided to investigate whether similar periodic behav- 
ior could also be found in the photometric data at our disposal. 
A correction of the photometry for the eclipse light curve can- 
not be carried out reliably since there is no accurate model light 
curve so far Using an improper model could result in the de- 
tection of spurious periodicities. Therefore we only used the 
photometric observations only from the eclipse period after the 
second photometric contact. The third contact had not yet oc- 
curred at the time of writing. We further assumed that the only 
brightness variations of e Aur in the UBV passbands were due 
to the intrinsi c variability of the only partially eclipsed F star 
(Kloppenbor g et alj|2010l) . i.e., we assumed the eclipse imposed 
no further photometric variability during this period. With these 
assumptions, we only needed to shift the two datasets of UBV 
photometry to a common zero point. We computed average val- 
ues for both datasets in all three passbands and subtracted them 
from the observed values. Analysis of these data by PERJOD04 
revealed that the observations in all filters can also be described 
by two harmonic functions with both periods nearly identical for 
all filters. The main period was again ~66 days. 

There are small differences in the main and secondary peri- 
ods found by this analysis, but these can be interpreted as a con- 
sequence of fitting data with finite noise on a finite time interval. 
Considering that the differences in all periods imply phase dif- 
ferences of only a few percent of a cycle at worst, we simply 
averaged the individual periods and their respective phase shifts. 
These mean quantities were subsequently held fixedm, while the 
fit to the observational data was repeated in order to determine 
their amplitudes. The resulting periods and phase shifts (with re- 
spect to the phase of the main RV period) are summarized in 
Tableg] 



12 



p. Chadima et al.: Spectral and photometric analysis of e Aurigae prior to and during the 2009-2011 eclipse 



The main period was found to be 66.21 days, and it is com- 
mon to all analyzed quantities, which proves that it is real. There 
is a phase shift of 0.808 between the RVs and the CIs and 0.323 
between the RVs and the UBV photometry. The secondary peri- 
ods differ among the various observables, but as already men- 
tioned, related observables share common secondary periods, 
which also suggests these periods are probably real. Figure [14] 
illustrates the final period fit for the RVs of the Fe II lines, for 
the CIs of the 6416.919 A Fe II line, and for the B photometric 
passband. 



periodicity. Attempts to extrapolate this periodic solution to 
the out-of-eclipse observations failed. The main period of the 
eclipse variability was found to be the same (66.21 days) for 
the radial velocity, central intensity, and photometric vari- 
ations. However, these observables did not vary in phase, 
but exhibited different phase shifts among the variables. The 
lack of a periodic solution to the out-of-eclipse data appears 
to rule out regular stellar pulsations as the cause of this vari- 
ability. The interpretation of the out-of-eclipse variability re- 
mains uncertain. 



6. Conclusions 

Our spectroscopic and photometric analysis of the e Aur led to 
the following conclusions as listed below. Since this study does 
not include observations from the end of the current eclipse and 
from post-eclipse phases, we intend to update this analysis once 
these phases have been covered by data from our ongoing obser- 
vations. 

1 . We attempted to disentangle the e Aur spectra and recover 
the spectra of both components, taking the unknown mass 
ratio q - Ml /M2 as a free parameter, but this effort was un- 
successful. We conclude that there is no secondary spectrum 
hidden in the prominent primary spectrum. However, it is 
possible that a secondary spectrum could remain hidden if 
the e Aur companion was itself a binary system. 

2. We derived a mean out-of-eclipse spectrum of the primary 
and compared it with a grid of synthetic spectra for 'classi- 
cal' F stars, using plane-parallel, LTE ATLAS model atmo- 
spheres. These analyses indicated that the F star is a low- 
gravity object, but a good spectral match was not found. 
However, this analysis should be repeated using more ap- 
propriate NLTE, spherically symmetric model atmospheres. 

3. We interpreted the complex behavior of the Ha line dur- 
ing eclipse as arising from additional absorption in an ex- 
tended 'atmosphere' around the companion's disk. We sub- 
tracted the mean out-of-eclipse Ha line profile from the pro- 
files observed during eclipse and obtained symmetric two- 
component absorption profiles. The behavior of these line 
profiles during the course of the eclipse agreed qualitatively 
with w hat is expected from the model presented by Huang 
(Il965h . This observational result, when combined with mod- 
eling of synthetic spectra for various eclipse phases, should 
serve to constrain the density and velocity structure of the 
companion's disk. 

4. A systematic decrease in the central intensity of the Ha line 
('spectroscopic eclipse' ) began ~3 years before photomet- 
ric first contact, owing additional circumstellar absorption. 
At the onset of spectroscopic eclipse, the line of sight to the 
F star lies far (at an angle of ~85°) from the mideclipse sight- 
line. This implies that circumstellar material associated with 
the secondary must extend well beyond the occulting disk re- 
sponsible for the photometric eclipses. In contrast, the sight- 
line at the start of photometric eclipse (first contact) lies only 
~20° away from the mideclipse line of sight. 

5. We searched for periodic behavior in the variable, out-of- 
eclipse Ha line, but did not find any periodicity. This im- 
plies that the distribution of circumstellar material around 
the primary (which may be a source of the Ha emission) is 
complex. 

6. During eclipse, both the Si II and the Fe II lines we exam- 
ined showed periodic variations in their radial velocities and 
central intensities. The UBV photometry showed a similar 
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